Magnetic field controlled electron transport in a thin cylinder 
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Abstract 

We explore electron transport in a thin cylinder, attached to two semi-infinite one-dimensional metallic 
electrodes, in the presence of both longitudinal and transverse magnetic fluxes. A simple tight-binding 
model is used to describe the system, where all the calculations are performed in the Green's function 
formalism. Quite surprisingly it is observed that, typical current amplitude oscillates as a function 
of the transverse magnetic flux, associated with conductance-energy characteristics, showing N(f)Q flux- 
quantum periodicity, where N and (j)Q (— ch/e) correspond to the system size and elementary flux- 
quantum respectively. The analysis might be helpful in fabricating mesoscopic switching devices, where 
a particular response can be delayed by tuning the system size N. 
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1 Introduction 

The sensitivity of electronic transport in low- 
dimensional systems like quantum dots [HE], meso- 
scopic loops [H H], quantum wires O [6], etc., on 
their geometry makes them truly unique in offering 
the possibility of electron transport in a very tun- 
able way. These systems provide several anomalous 
characteristics in electron transport due to their 
reduced system dimensionality and quantum con- 
finement. In the present age of nanoscience and 
technology, electronic devices made from quantum 
confined model systems are used extensively in fab- 
rication of nano-electronic circuits where electron 
transport becomes predominantly coherent [7l [8]. 
The progress of theoretical description in a bridge 
system has been followed based on the pioneer- 
ing work of Aviram and Ratner Later, the 
actual mechanisms underlying such transport be- 
come much more clearly resolved after several nice 
experimental observations [lOl [Til [12] in different 
bridge systems. Though in literature many the- 
oretical as well as experimental papers on elec- 
tron transport are available, yet lot of controver- 
sies are still present between the theory and ex- 
periment, and the complete knowledge of the con- 
duction mechanism in this scale is not very well 
established even today. Several key factors are 
there which control the electron transport across a 
bridge system, and all these effects have to be taken 
into account properly to characterize the trans- 
port. For our illustrative purposes, here we de- 
scribe very briefiy some of these factors. (I) The 
quantum interference effect [H [H [TH [THl [H] of 
electron waves passing through different arms of a 
conducting element which bridges two semi-infinite 
one-dimensional metallic electrodes, viz, source and 
drain becomes the most significant issue in electron 
transport in a bridge system. (II) The coupling of 
the electrodes with the bridging material provides 
an important signature in the determination of cur- 
rent amplitude 13J. The understanding of this cou- 
pling under non-equilibrium condition is a major 
challenge in this particular field, and we should take 
care about it in fabrication of any electronic device. 
(Ill) The geometry of the conducting material be- 
tween the two electrodes itself is an important issue 
to control the electron transmission. To empha- 
size it, Ernzerhof et al. [18j have predicted several 
model calculations and provided some new signif- 
icant results. (IV) The dynamical fluctuation in 
the small-scale devices is another important factor 
which plays an active role and can be manifested 



through the measurement of shot noise, a direct 
consequence of the quantization of charge [TU [20] . 
Beside these factors, several other parameters of the 
Hamiltonian that describe a system also provide sig- 
nificant effects in the determination of the current 
across a bridge system. 

In the present paper, we will investigate the 
electron transport properties of a thin cylinder 
(see Fig. [T]), attached to two semi-infinite one- 
dimensional metallic electrodes, in the presence of 
both longitudinal and transverse magnetic fluxes, 0; 
and 0t respectively. A simple tight-binding model 
is used to illustrate the system, where all the cal- 
culations are performed in the Green's function for- 
malism. Our numerical study shows that, typical 
current amplitude across the cylinder oscillates as 
a function of the transverse magnetic flux 0(, as- 
sociated with the radial magnetic field Br, show- 
ing N(f)Q flux-quantum periodicity instead of sim- 
ple 00-periodicity, where N and (fio {— ch/e) corre- 
spond to the size of the cylinder and the elementary 
flux-quantum respectively. This oscillatory behav- 
ior provides an important signature in this particu- 
lar study. It is observed that, for a fixed bias voltage 
V applied by the electrodes, we can achieve a simi- 
lar response for different cylinders by controlling the 
transverse magnetic flux associated with Br- Thus, 
for the two cylinders of different sizes, similar re- 
sponse will be obtained for two different values of 
Br, depending on the system size of the cylinder. 
For smaller cylinder, the response is obtained for 
lower value of Br, while this similar response will be 
obtained for some higher value of Br for the larger 
cylinder. This clearly manifests that a particular 
response can be delayed by tuning the system size 
N quite significantly. This aspect may be utilized 
in designing a tailor made switching device, and to 
the best of our knowledge, this phenomenon has not 
been addressed earlier in the literature. 

We organize the paper as follow. Following the 
introduction (Section 1), in Section 2, we present 
the model and the theoretical formulations for our 
calculations. Section 3 discusses the significant re- 
sults, and finally, we summarize our results in Sec- 
tion 4. 

2 Model and the theoretical 
description 

This section follows the methodology for the cal- 
culation of the transmission probability (T), con- 
ductance (g) and current (/) in a small cylinder by 
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using the Green's function technique. 

We begin by referring to Fig. [U where a cyHn- 
der is attached to two semi-infinite one-dimensional 
metaUic electrodes, viz, source and drain. The 
cylinder is subjected to a longitudinal magnetic flux 

, and a transverse magnetic flux (j)t . For our illus- 




Figure f: (Color online). Schematic view of a 
mesoscopic cylinder attached to two semi-infinite 
one-dimensional metallic electrodes and subjected 
to both the longitudinal and transverse magnetic 
fluxes and (j)t respectively. Filled circles corre- 
spond to the position of the atomic sites. 

trative purpose here we consider the simplest pos- 
sible cylinder, where only two isolated one-channel 
rings are connected by some vertical bonds. The 
transverse magnetic flux (pt is expressed in terms 
of the radial magnetic field Br by the relation 
4>t = BrLd, where the symbols L and d correspond 
to the circumference of the ring and the hight of the 
cylinder respectively. 

At sufficient low temperature and bias voltage, 
we use the Landauer conductance formula [2TJ [22] 
to calculate the conductance g of the cylinder which 
can be expressed as, 



where T becomes the transmission probability of an 
electron through the cylinder. It can be expressed 
in terms of the Green's function of the cylinder and 
its coupling to the two electrodes by the relation [211 

T^Tr [TsG'^dG''^] (2) 

where and G° are respectively the retarded and 
advanced Green's functions of the cylinder includ- 
ing the effects of the electrodes. The parameters 
Fs and To describe the coupling of the cylinder 
to the source and drain respectively, and they can 



be defined in terms of their self-energies. For the 
complete system i.e., the cylinder with the two elec- 
trodes the Green's function is defined as, 

G = {e-H)-' (3) 

where e — E + irj. E is the injecting energy of the 
source electron and 77 gives an infinitesimal imagi- 
nary part to e. Evaluation of this Green's function 
requires the inversion of an infinite matrix as the 
system consists of the finite cylinder and the two 
semi-infinite electrodes. However, the entire system 
can be partitioned into sub-matrices corresponding 
to the individual sub-systems and the Green's func- 
tion for the cylinder can be effectively written as, 

Ge = (e - i?c - - ^dV' (4) 

where He is the Hamiltonian of the cylinder which 
can be written in the tight-binding model within 
the non-interacting picture like, 

N N 

1=1 1=1 

<ij> 

<ij> 
N 

i=l 

In the above Hamiltonian {He), ef's (e[^'s) are the 
site energies in the lower (upper) ring, c^^ (c[^^) 
is the creation operator of an electron at site i 
in the lower (upper) ring, and the corresponding 
annihilation operator for this site i is denoted by 
(cf). The symbol vf" {v^ ) gives the nearest- 
neighbor hopping integral in the lower (upper) ring, 
while the parameter vt corresponds to the trans- 
verse hopping strength between the two rings of 
the cylinder. 9i and 62 are the two phase fac- 
tors those are related to the longitudinal and trans- 
verse fluxes by the expressions, 9i — 27T<j)i/NipQ and 
02 — 7r(/)(/7V0o, where N represents the total num- 
ber of atomic sites in each ring and 0o = ch/e 
is the elementary flux-quantum. Similar kind of 
tight-binding Hamiltonian is also used to describe 
the two semi-infinite one-dimensional perfect elec- 
trodes where the Hamiltonian is parametrized by 
constant on-site potential eo and nearest-neighbor 
hopping integral to- In Eq. [31 E5 — hg^gshsc and 
'^D — hucgoh^Dc are the self-energy operators due 
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to the two electrodes, where gs and go correspond 
to the Green's functions of the source and drain 
respectively, hsc and huc are the coupling matri- 
ces and they will be non-zero only for the adjacent 
points of the cylinder, and the electrodes respec- 
tively. The matrices and To can be calculated 
through the expression. 



S{D) 



(6) 



where S^^^^ and S^^^^ are the retarded and ad- 
vanced self-energies respectively, and they are con- 
jugate with each other. Datta et. al. |^ have 
shown that the self-energies can be expressed like 
as, 



^S{D) 



Trn. 



= ^s{D) - i^s{D) (7) 
where h^s(D) are the real parts of the self-energies 
which correspond to the shift of the energy eigenval- 
ues of the cylinder and the imaginary parts l^s(D) of 
the self-energies represent the broadening of these 
energy levels. This broadening is much larger than 
the thermal broadening and this is why we restrict 
our all calculations only at absolute zero temper- 
ature. All the informations about the cylinder-to- 
electrode coupling are included into these two self- 
energies. 

The current passing across the cylinder is de- 
picted as a single-electron scattering process be- 
tween the two reservoirs of charge carriers. The 
current / can be computed as a function of the ap- 
plied bias voltage V through the relation ^21J, 

=-+eV/2 

T{E,V)dE (8) 

-eV/2 

where Ep is the equilibrium Fermi energy. For 
the sake of simplicity, we assume that the entire 
voltage is dropped across the cylinder-electrode in- 
terfaces and this assumption doesn't greatly affect 
the qualitative aspects of the I-V characteristics. 
Such an assumption is based on the fact that, the 
electric field inside the cylinder especially for small 
cylinders seems to have a minimal effect on the 
conductance-voltage characteristics. On the other 
hand, for quite larger cylinders and high bias volt- 
ages the electric field inside the cylinder may play 
a more significant role depending on the internal 
structure and size of the cylinder [25 , yet the effect 
is quite small. 

3 Results and discussion 

To reveal the basic mechanisms of and 0f on the 
electron transport, here we present all the results 



only for the non-interacting electron picture. With 
this assumption, the model becomes quite simple 
and all the basic features can be well understood. 
Another realistic assumption is that, we focus on 
the perfect cylinders only i.e., the site energies are 



taken as ef 



for all i. As illustrative pur- 



poses, we parametrize all the hopping integrals (w; , 
vY and vt) of the cylinder by a single parameter v, 
viz, we set vf" = vY = vt = v. In the present arti- 
cle, we concentrate our study of electron transport 
in the two separate regimes. These are so-called 
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Figure 2: (Color online). g-E characteristics for a 
thin cylinder with N = 12, (pi = 0.2 and <j)t = 2. (a) 
weak-coupling limit and (b) strong-coupling limit. 

the weak-coupling and strong-coupling regimes re- 
spectively. The weak-coupling regime is specified by 
the condition t^s,d} << On the other hand, the 
condition T{5,_d} ^ v is used to denote the strong- 
coupling regime. Here, the parameters ts and to 
correspond to the hopping strengths of the cylinder 
to the source and drain respectively. In these two 
limiting cases, we choose the values of the differ- 
ent parameters as follow: ts = td = 0.5, v = 2.5 
(weak-coupling) and ts = td = 2, v = 2.5 (strong- 
coupling). The on-site energy eo is set to (we can 
take any constant value of it instead of zero, since 
it provides only the reference energy level) for the 
electrodes, and the hopping strength to is taken as 3 
in the two semi-infinite metallic electrodes. For the 
sake of simplicity, we set the Fermi energy Ep — 
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and choose the units where c = e — h = I. 

Let us begin our discussion with the variation 
of the conductance g as a function of the inject- 
ing electron energy E. As iUustrative examples, in 
Fig. [2] we plot the conductance-energy (g-E) char- 
acteristics for a thin cylinder with N = 12, consid- 
ering the parameters (f>i = 0.2 and (j)t — 2. Fig- 
ures mja) and (b) correspond to the results for the 
weak- and strong-coupling regimes respectively. It 
shows that, in the limit of weak-coupling, the con- 
ductance exhibits fine resonance peaks for some par- 
ticular energies, while for all other energies the con- 
ductance almost disappears. Out of these resonance 
peaks, only for the two energies the conductance ap- 
proaches the value 2. Thus, for these energies the 
transmission probability T becomes unity, since we 
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fects, and hence the cylinder to electrodes interface 
structure. Now all these resonance peaks are as- 
sociated with the energy eigenvalues of the cylin- 
der, and accordingly, we can say that the conduc- 
tance spectrum manifests itself the electronic struc- 
ture of the cylinder. In the resonance spectrum we 
get more peaks in the presence of both (pi and (j)t 
compared to those as obtained when both these two 
fluxes are identically zero. This is due to the fact 
that the fluxes and 4>t remove all the degenera- 
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Figure 3: (Color online). Typical conductance as a 
function of for a thin cylinder with N = 12 in 
the strong-coupling limit, where the parameter (j)t 
is set to 2. The conductance is calculated at the 
energy E = 2.1. 

get the relation g = 2T from the Landauer con- 
ductance formula (see Eq. [T] with e — h = 1). 
On the other hand, for other resonances, the con- 
ductance does not reach to 2 any more, and gets 
much reduced value. This is due to the quantum 
interference effect of the electronic waves passing 
through the different arms of the cylinder, and it 
can be interpreted as follow. During the motion 
of the electrons from the source to drain through 
the cylinder, the electron waves propagating along 
the different possible pathways can get a phase shift 
among themselves, according to the result of quan- 
tum interference. Therefore, the probability ampli- 
tude of getting an electron across the cylinder either 
becomes strengthened or weakened. This causes 
the transmittance cancellations and provides anti- 
resonances in the conductance spectrum. Thus it 
can be emphasized that the electron transmission 
is strongly affected by the quantum interference ef- 
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Figure 4: (Color online) . Typical conductances as a 
function of (pt for three different cylinders with (a) 
iV = 4, (b) iV = 6 and (c) iV 8 respectively, in 
the limit of weak-coupling. The energy E is fixed 
at 2.1 and the parameter 0/ is set to 0.2. 

cies in the energy levels which provide more reso- 
nance peaks in the conductance spectrum. With 
these features, we get additional one feature when 
the coupling strength of the cylinder to the elec- 
trodes increases from the low regime to high one. 
In the limit of strong cylinder-to-electrode coupling, 
all these resonances get substantial widths com- 
pared to the weak-coupling limit. The results are 
shown in Fig.[2I^b). The contribution for the broad- 
ening of the resonance peaks in this strong-coupling 
limit appears from the imaginary parts of the self- 
energies Eg and E^i respectively [21], as mentioned 
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earlier. Thus by tuning the couphng strength, we 
can get the electron transmission across the cylin- 
der for wider range of energies and it provides an 
important signature in the study of current-voltage 
characteristics. 

In order to express the dependence of both the 
two different fluxes on the electron transport much 
more clearly, in Figs. [3] and H] we display the varia- 
tion of the typical conductances as a function of (f>i 
and (pt respectively. Figure [3] shows the variation of 
the typical conductance (gtyp) with the longitudinal 
magnetic flux (f>i in the limit of strong-coupling for 
a cylinder considering = 12. Here we calculate 
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Figure 5: (Color online). I-V characteristics for a 
thin cylinder with N = 6, (pi — 0.2 and (pt = 2. (a) 
weak-coupling limit and (b) strong-coupling limit. 

gtyp for the energy E = 2.1, when the parameter 
ipt is set to 2. It shows that the typical conduc- 
tance exhibits oscillatory behavior with (pi showing 
00 (= 1, since c = e = /i = lin our present formula- 
tion) flux-quantum periodicity. This (po periodicity 
is well known in the literature. The remarkable be- 
havior is observed only for the case when we apply 
the transverse magnetic flux in the cylinder. For 
the illustrative purposes, in Fig. 2] we plot the varia- 
tion of the typical conductance of the three different 
cylinders in the limit of weak-coupling, where (a), 
(b) and (c) correspond to the results for the cylin- 
ders with = 4, 6 and 8 respectively. In this case, 
we set the energy E = 2.1 and fix <pi to 0.2. Quite 



interestingly it is observed that, the typical conduc- 
tance varies periodically with the flux (pt showing 
N<po flux-quantum periodicity, instead of the tra- 
ditional 00 flux-quantum periodicity as observed in 
the previous case. The detailed investigation re- 
veals that the nature of the gtyp vs. (pt characteris- 
tic within a single period of any one of these three 
cylinders is quite similar to that of the other two 
cylinders within their respective periods. Thus the 
positions of the resonance peaks shift appropriately 
with the increase of the system size N to make the 
nature invariant. This is really a very interesting 
phenomenon and can be utilized to get a delayed 
response just by tuning the system size N. 

All these features of electron transport become 
much more clearly visible by studying the current- 
voltage (I-V) characteristics. The current across 
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Figure 6: (Color onhne) . Typical current amplitude 
as a function of (pi for a thin cylinder with TV = 6 
in the strong-coupling limit, where the parameter 
(pt is set to 2. The current is calculated at the bias 
voltage V = 1.95. 

the cylinder is computed from the integration pro- 
cedure of the transmission function T as prescribed 
in Eq. [H The transmission function varies exactly 
similar to that of the conductance spectrum, dif- 
fer only in magnitude by the factor 2 since the re- 
lation g = 2T holds from the Landauer conduc- 
tance formula Eq. [TJ As representative examples, 
in Fig. [5] we plot the current-voltage characteristics 
for a thin cylinder with N — 6, (pi — 0.2 and cpt = 2, 
where (a) and (b) represent the results for the weak 
and strong cylinder-to-electrode coupling limits re- 
spectively. In the limit of weak-coupling, the cur- 
rent exhibits staircase-like structure with fine steps 
as a function of the applied bias voltage. This is 
due to the existence of the sharp resonance peaks 
in the conductance spectrum in this limit of cou- 
pling, since the current is computed by the integra- 
tion method of the transmission function T. With 
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the increase of the apphed bias voltage, the elec- 
trochemical potentials on the electrodes are shifted 
gradually, and finally cross one of the quantized 
energy levels of the cylinder. Accordingly, a cur- 
rent channel is opened up and the current-voltage 
characteristic curve provides a jump. On the other 
hand, for the strong cylinder-to-electrode coupling, 
the current varies almost continuously with the ap- 
plied bias voltage and achieves much larger ampli- 
tude than the weak-coupling case. The reason is 
that, in the limit of strong-coupling all the energy 
levels get broadened which provide larger current in 




Figure 7: (Color online). Typical current ampli- 
tudes as a function of 4>t for two different cylinders 
with (a) = 4 and (b) = 6 respectively, in the 
limit of strong-coupling. The voltage V is fixed at 
0.7 and the parameter 0; is set to 0.3. 



voltage y = 1.95. It predicts that the current am- 
plitude shows an oscillatory behavior with i^;, giv- 
ing (pQ flux-quantum periodicity. This is completely 
analogous to the variation of gtyp with (pi as de- 
scribed earlier. The interesting feature is observed 
only when we study the variation of Ityp with 4>t, 
instead of 0/ . The results are shown in Fig. [7l where 
(a) and (b) correspond to the variation of the cur- 
rent amplitudes for the cylinders with A^ = 4 and 
6 respectively. The typical current amplitudes are 
calculated in the limit of strong-coupling for the 
bias voltage V = 0.7, when the parameter (f>i is 
set to 0.3. Quite interestingly we see that, the cur- 
rent amplitude varies periodically with (f>t providing 
N(j)Q flux-quantum periodicity, instead of the con- 
ventional periodicity. Thus for the cylinder with 
A^ = 4, the current shows 40o (= 4) periodicity, 
while it becomes 60o (= 6) for the cylinder with 
A^ = 6. From these results we clearly observe that 
the variation of Ityp with (pt within a period for the 
cylinder with A^ = 6 is quite similar to that of the 
cylinder with A^ = 4 within its single period. This is 
just the replica of the gtyp versus 4>t characteristics 
which we have described earlier. Thus for a fixed 
bias voltage, we can achieve the similar response 
for different values of A just by tuning the flux (pt 
associated with the radial magnetic field Br- This 
phenomenon can be utilized as follow. Let us con- 
sider two such cylinders of different sizes. For these 
two cylinders, a particular response is obtained for 
the two different values of B^, depending on the 
size of these cylinders. Therefore, for the cylinder 
of smaller size the similar response is achieved at 
the lower value of Br compared to the cylinder of 
larger size. This clearly manifests that a particu- 
lar response can be delayed by tuning the size A of 
the cylinder quite significantly. This aspect is re- 
ally a very interesting, and provides a signature for 
manufacturing mesoscopic switching devices. 



the integration procedure of the transmission func- 
tion T . Thus by tuning the strength of the cylinder- 
to-electrode coupling, we can achieve very large cur- 
rent amplitude from the very low one for the same 
bias voltage V . 

Finally, we concentrate our study on the varia- 
tion of the typical current amplitude {hyp) with the 
magnetic fluxes pi and pt to illuminate the signifi- 
cant effects of these fluxes on the electron transport. 
In Fig. [51 we display the variation of Ityp as a func- 
tion of 4>i for a thin cylinder with A^ = 6 in the limit 
of strong-coupling, where the parameter (f>t is fixed 
to 2. The typical current is computed at the bias 



4 Concluding remarks 

In conclusion, we have addressed the electron trans- 
port in a thin cylinder, attached to two semi-infinite 
one-dimensional metallic electrodes, in the pres- 
ence of both longitudinal and transverse magnetic 
fluxes. We have used a simple tight-binding model 
to describe the system where all the calculations 
have been done in the Green's function formalism. 
Quite interestingly we have observed that the typ- 
ical current amplitude oscillates as a function of 
the transverse magnetic flux pt, associated with 
the conductance-energy characteristics, providing 
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N(I)Q flux-quantum periodicity. This feature is com- 
pletely different from the traditional oscillatory be- 
havior like as we have got for the case of Ityp versus 
(pi characteristic. Our results have predicted that 
a particular response can be delayed by tuning the 
size of the cylinder N. This aspect may be utilized 
in designing a tailor made mesoscopic switching de- 
vice. 

This is our first step to describe how the electron 
transport in a thin cylinder can be controlled very 
nicely by means of the longitudinal and transverse 
magnetic fluxes. We have made several realistic as- 
sumptions by ignoring the effects of the electron- 
electron correlation, disorder, temperature, finite 
width of the electrodes, etc. Here we discuss very 
briefly about these approximations. The inclusion 
of the electron-electron correlation in the present 
model is a major challenge to us, since over the last 
few years many people have studied a lot to incor- 
porate this effect, but no such proper theory has yet 
been developed. In this work, we have presented all 
the results only for the ordered systems. But in real 
samples, the presence of impurities will affect the 
electronic structure and hence the transport proper- 
ties. The effect of the temperature has already been 
pointed out earlier, and, it has been examined that 
the presented results will not change significantly 
even at finite temperature, since the broadening of 
the energy levels of the cylinder due to its coupling 
with the electrodes will be much larger than that of 
the thermal broadening |21j . The other important 
assumption is that here we have chosen the linear 
chains instead of wider leads, since we are mainly 
interested about the basic physics of the cylinder. 
Though the results presented here change with the 
increase of the thickness of the leads, but all the 
basic features remain quite invariant. Finally, we 
would like to say that we need further study in this 
system by incorporating all these effects. 
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